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Abstract 

We introduce an exactly solvable example of timelike geodesic motion and geodesic deviation 
in the background geometry of a well-known two-dimensional black hole spacetime. The effective 
potential for geodesic motion turns out to be either a harmonic oscillator or an inverted harmonic 
oscillator or a linear function of the spatial variable, corresponding to the three different domains 
of a constant of the motion. The geodesic deviation equation also is exactly solvable. The cor- 
responding deviation vector is obtained and the nature of the deviation is briefly discussed by 
highlighting a specific case. 
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I. INTRODUCTION 



Learning to solve the geodesic and geodesic deviation equations is an integral part of a 
first course on general relativity. However, most texts on general relativity^ do not contain a 
sufficient number of solvable examples that illustrate the behavior of geodesies in a nontrivial 
geometry. This lack is largely due to the fact that the geodesic equations as well as the 
geodesic deviation equations are not easy to solve exactly in most of the standard geometries 
of interest. A traditional example is the case of the line element on the surface of a two- 
dimensional sphere for which the geodesic and the geodesic deviation equations can be 
solved analytically. Some other exactly solvable examples in two dimensions include Rindler 
spacetime, the line element on a cylinder and the hyperbolic plane. 

In this article, we will discuss a two-dimensional example (with a diagonal line element 
of Lorentzian signature ( — h), i.e. of the form ds 2 = g^dt 2 + gudx 2 , with g 0Q everywhere 
negative and gn, positive) , for which, like the sphere, there are exact solutions for geodesies 
as well as geodesic deviation. The line element for this example arises in the context of the 
stringy black holes.— In fact, this is the geometry of the first stringy black hole discovered in 
the early nineties by Mandal et al.— We will not go into the details of how this line element 
is obtained. Interested readers may read the original references^ to enlighten themselves on 
this aspect. For our purposes here, the line element represents a geometry that is given to 
us, and we see that there are exact solutions for geodesies and geodesic deviation for test 
particles moving in this geometry. 



II. THE LINE ELEMENT 



We begin by writing down the line element of interest 4 : 

, ? Ms , 9 kdr 2 

ds 2 = -( l ~-) dt + 4r2 (i-M) - C 1 ) 

where M is equivalent to GM/c 2 with the choice of units G = 1 and c = 1. Similarly k is 
required to have dimensions of length squared in order to render the metric coefficient g rr 
dimensionless. 

The domain for the radial coordinate is r > M and t is allowed to vary from — oo < t < oo. 
The physical meaning of the constants M and k can be explained in the context of two- 
dimensional string theory. For completeness (and also with the risk of introducing some 



unexplained jargon), we mention that M is related to mass and k to the central charge 
parameter. 

Several purely geometrical facts about this line element are worth noting. At r = M, 
g u — > and g rr —>■ oo. If r < M, the signs of g u and g rr are interchanged. Readers familiar 
with some general relativity (say, the Schwarzschild solution) would call r = M the horizon 
and thereby, characterize the line element as a black hole provided there exists a singularity 
(which is the case here, as we shall see later) inside the horizon (that is, in the region r < M). 
As the title of this article suggests, the geometry represents a two-dimensional black hole. As 
is common knowledge, a black hole is a region from where nothing can escape. If an object 
falls into a black hole (that is, the object crosses the horizon), the object ends up getting 
stretched and torn apart by the gravitational forces within. This distortion is quantified 
through the deviation of geodesies. The ultimate fate is collapse into a singularity where 
the matter density is infinitely large, thereby being a region of extreme curvature. We shall 
concern ourselves with the r > M region. 

The pathology in the line element (that is, g rr —>■ oo ) at r = M is entirely a coordinate 
artifact (a poor choice of coordinates). In addition, it is worth noting, that unlike the case for 
similar but four- dimensional line elements where the horizon is a sphere or some other two- 
dimensional surface, here it is a point and falling into a black hole is manifest in crossing this 
point and entering the region of no return. The spacetime here is undoubtedly very similar to 
the two-dimensional r-t section of the static, spherically symmetric four-dimensional black 
holes (like the Schwarzschild metric). Many of the features of black holes can be understood 
via this example. In fact an extensive research literature exists on two-dimensional black 
holes where interesting phenomena such as Hawking radiation and information loss have 
been discussed (see Ref. 3 and references therein and Ref. O). 

The line element in Eq. (^) also can be written in the following form 3 : 

ds 2 = -tanh. 2 ( ^=)dt 2 + dx 2 , (2) 

where we have used the coordinate transformation x = ^cosh _1 (|j — 1) > 0. Notice that 
the singular nature of the metric function gu at r = M is no longer there in this form of 
the line element. A third form of the line element is the Kruskal extension or the maximal 
analytical extension done exactly in the same way as for the Schwarzschild metric. (For 
details on how to construct Kruskal coordinates see any standard text on general relativity.) 



If we introduce the coordinates u and v defined by: 

\fk 

u = t- — ln|r - M\ (3a) 

v = t+ ^-ln|r - M\, (3b) 

and U = —y/kexp(—u/\/~k) J V = \/~kexp(v /Vk), we can rewrite the line element as: 

2 dUdV 
= ~ m uv- v*J 

Vk k 

In Eq. (J3J) the coordinates u and v run from — oo to oo and — oo < U < 0, < V < oo. It 
is possible to further rewrite the line element using coordinates 2T = V + U, 2X = V — U . 
The Kruskal form of the line element is an 'extension' over the whole domain of r and t The 
geometry of this Kruskal extension has features similar to that of the Schwarzschild. 



III. TIMELIKE GEODESICS 

Our main goal is to find timelike geodesies (that is, trajectories of test particles of nonzero 
rest mass) in the line element given in Eq. (0). The equations and constraint for timelike 
geodesies are given as: 

x^ + T» x x u x x = (5a) 
g^afx" = -1, (5b) 

where are the coordinates (here r and t), and \i takes on values and 1 for t and r 
respectively. Equation (5b) is the timelike constraint on the trajectories 

The differentiation (indicated by dots) in Eq. (0) is with respect to the affine parameter 
(proper time) A. The Christofel symbol T^ x is given as: 

Kx = \9 w \9 P v,\ + 9 P x,u - guxj, (6) 

where a comma denotes ordinary differentiation with respect to coordinate x^.The nonzero 
components of T^ x for the line element in Eq. (Q) are 

r ?o = r oi = 2r(r _ M) ( 7a ) 

r^o = (Tb) 
r 1 - {2r ~ M) (7c) 

11 ~ 2r(r-M) [<C) 



The reader might ask why aren't we interested in null geodesies (that is, trajectories of 
test particles having zero rest mass, such as light rays) ? In two dimensions null geodesies are 
the same as those for the flat Minkowski line element. Note that any two-dimensional line el- 
ement (with Euclidean or Lorentzian signature) can be written as ds 2 = Q 2 (x,t)[— dt 2 + dx 2 ], 
where Q 2 (x, t) is a nonzero, positive definite function of x and t and is called the conformal 
factor. Null geodesies require that g^ v x^x u = 0, which implies that in two dimensions, null 
geodesies for Minkowski spacetime are the same as those for any other nontrivially curved 
spacetime. 

The timelike constraint in Eq. (JHJ) for the line element of interest here (Eq. (0)), is given 

as: 

_(1 _ E)i 2 + j^-^r 2 = -1. (8) 

V r ; 4r 2 (l - f ) V ; 

The geodesic equations (one each for r and t) are 

* = 2(TTM) (9) 

Ir2r-M-|. 2 E 2 Mr 
r_ 2l 2 -M/ + 2k(r — M) ' ( } 

where in Eq. (JSJ) we have substituted for t from Eq. (JSJ). E in the above is a constant of the 

motion and arises due to the absence of the coordinate t in the metric coefficients. If we 

substitute t from Eq. Q in the timelike constraint in Eq. (|5|). we obtain 

r 2 = h(E 2 - A)r 2 + 4Mr] = -V{r) + E , (11) 

K 

where Eq is a constant and V(r) is like an effective potential given by: 

V(r) = E - h(E 2 - 4)r 2 + 4Mr}. (12) 

rC 

Notice that the effective potential is like that of the harmonic oscillator for E 2 < 4 and 
like that of an inverted harmonic oscillator for E 2 > 4. For E 2 = 4 it is linear with a 
negative slope. We can make a coordinate change (translation) in order to see the harmonic 
oscillator /inverted harmonic oscillator forms explicitly. The potentials for the three cases 
are shown in Fig. 1. 

If we substitute r and its derivative with respect to A from the timelike constraint (jllj) 
in Eq. (fTUj) for r, the equation is satisfied identically. Therefore, it is enough to just obtain 



r(A) by integrating the timelike constraint for the three cases (E 2 < 4, E 2 = 4, E 2 > 4) and 
obtain r(A) and subsequently t(A). These results are given below. 
Case 1: E 2 < 4. 

r(A) = — (1 + cosaA) (13a) 

£ r , 1 , o + tancA-, . 

'W^V^F^ 1 ' < 13b) 
where q 2 = (1 + /)/(l - /), A = 4 - £ 2 , a = ^/A/k, c = a/2, f = 1 - A/2, and B = AM. 
Case 2: £ 2 = 4. 

r(A) = ^A 2 (14a) 

t (A) = A + ^(^). (14b) 

Case 3: £ 2 > 4. 

r(A) = — (-1 + coshaA) (15a) 
E Vs o, tanhcA — o-, . 

'W^ + ^I^aT^ (15b) 

where q 2 = (f - l)/(/ + 1), A = E 2 - 4, a = c = a/2, / = 1 + A/2, and 5 = 4M. 

In each case, the functional form of r(t) can be obtained by utilizing the parametric forms 
r(A), t(X). Figures 2-4 illustrate the above function. The origin of the coordinates in each 
of these figures is chosen to be r = 1 and t — 0. In particular, Fig. 4 shows the actual radial 
trajectory r(t) as a function of time. It is not possible to write r as an explicit function of 
t because of the nature of the functional dependence of t on A. However, we may write t(r) 
by inverting r(A) to A(r) and then using this function in t(A(r)). Notice that for E 2 < 4, 
the trajectory hovers around the horizon (extending only up to r = 2M = 2), reaching the 
value r — M — 1 only asymptotically. However, for E 2 > 4 there is no such restriction and 
the trajectory can extend from M to oo. 

It should be mentioned that the "discontinuities" in some of the plots arise because of 
the choice of the parameter (M = 1 implies that r > 1) and, consequently, the allowed 
domain of the independent variable. They are not real discontinuities of the functions, as is 
apparent from the solutions. 
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IV. GEODESIC DEVIATION 



We now consider the geodesic deviation equation in order to analyze the behavior of 
the spreading of the geodesic curves obtained in Sec. 11111 Before we discuss the geodesic 
deviation in the black hole line element, we briefly recall some of the basic notions of geodesic 
deviation. 5 

In a curved spacetime geodesies are not the straight lines of our Euclidean intuition. 
For instance, the geodesies on the surface of a sphere are the great circles. To understand 
how geodesies are spread out in a spacetime of nontrivial curvature, we consider a bundle 
of geodesies around a specific geodesic (called the central geodesic). At every point on 
this central geodesic we erect little normals. This normal is the connecting vector, so-called 
because as we move along the connecting vector, we reach neighboring geodesies. The reader 
might ask why normal deformations only? The answer is that tangential deformations are 
taken into account through the reparametrisation of the afnne parameter A, or, in other 
words, they do not deform the trajectory to a neighboring one. 

It is obvious that the nature of the connecting vector (or deviation vector) in a nontrivial 
spacetime depends crucially on the properties of the spacetime, or more importantly, as it 
turns out, on its curvature properties. For instance, on the surface of a sphere, if we look at 
two neighboring great circle arcs connecting the north and south poles, the distance between 
them increases, reaches a certain value at the equator, and then decreases untill it reaches 
a zero value at the other pole. The deviation vector essentially measures this effect, which 
is dictated by the properties of the Riemann tensor R^ Xp given by: 

DM _ pA* _ pM I pA* per _ p(U per 

Il u\p~ L up,\ L u\,p^ L a\ L vp L ap L uX- \ LU ) 

The Ricci tensor and the Ricci scalar are obtained from the above Riemann tensor by 
appropriate contractions with the metric tensor (Ricci tensor, R^ v = R* Xv an d Ricci scalar, 
R = g^Rpu)- 

We know from Einstein's general relativity that geometry is equivalent to the presence 
of a gravitational field. Thus the deviation of trajectories from each other is a measure of 
the relative gravitational force between the objects moving on these separate, neighboring, 
trajectories. We now state the geodesic deviation equation and solve it for the geodesies in 
the black hole line element, Eq. discussed here. 
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Given a geodesic curve, we identify the tangent and normal to it as the vectors and 
respectively. These satisfy the orthonormality conditions: 



g^ePe" = -1 (17a) 
g^n" = 1 (17b) 
g^n" = 0. (17c) 



We make the assumptions: 



e M = (t,r) (18a) 
n" = (f(r)r,g{r)i). (18b) 

Note that other choices can be made and thus the above definitions in Eq. (|18|) are by 
no means unique. If we substitute the assumptions in Eq. ()18j) into the orthonormality 
conditions Eq. (|17jl. we obtain the following relations for f(r) and g(r): 



f(r) = J J*- (19a) 

-goo 



9(r) = J=&. (19b) 
V 9u 

Therefore the normal vector n M takes the form: 

2r(l-^). , 
n 1 = V r U . (20b) 
v 

The equation for the normal deformations 77 (where 77 M = 7771^, 77^ being the deviation vector) 
is given by: 

■^ + R^ p ^n u e^ V = 0. (21) 

The nonzero components of Riemann tensor in the coordinate frame are given by: 

M 

-R0101 — -R1010 — — (22a) 
M 

-R0110 — -R1001 — ^~3- (22b) 

There is actually only one independent component of the Riemann tensor in two dimen- 
sions (recall that in n dimensions, the number of independent components is n 2 (n 2 — 1)/12). 
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The nonzero Ricci tensor components in the coordinate frame (R^ = R^ au ) and the Ricci 
scalar (R = g^R^) are 

_ 2M(r - M) 
K °° ~ yT2 

R " = M7^M) (23b) 

R = ^ (23c) 

The Riemann tensor components and Ricci scalar diverge as r — > 0. The geometry 
becomes singular as r — > 0. However, r — > M does not seem to be a preferred point (in 
comparison to r = say), a fact that re-emphasizes the point that the pathology at r = M 
in the line element in Eq. is largely a coordinate artifact. Of course, the special nature 
of r = M is that it is a horizon with the metric signature changing sign as we move from 
r > M to r < M. 

If we substitute the Riemann tensor components for the two-dimensional black hole line 
element and use the normal and tangent mentioned in Eq. , we find that the deviation 
equation reduces to the following equations for the three cases discussed earlier. 

Case 1: E 2 < 4. 

d 2 r] 



d\ 2 



- 2 sec 2 Xq = 0, (24) 



where A = ~\ 4 u^- A. The two linearly independent solutions to Eq. (|2l|) are 



Case 2: E 2 = 4. 



The solutions are 



Case 3: E 2 > 4. 



77(A) = CitanA + C 2 (AtanA + 1). (25) 
d 2 r] 2 

J ~ V = °- (26) 

,(A) = dA 2 + C 2 i (27) 



' 2cosech 2 A77 = 0, (28) 



dX 



2 



where A = ~W E2 fc 4 A. In this case, the solutions are 



77(A) = C 1 cothA + C 2 (AcothA-l) (29) 
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In Eqs. (|25|). ()27|). and ()29j) . Ci and C2 are separate arbitrary constants in each case. These 
constants should be determined by the initial conditions (that is, the value of the initial 
separation between a pair of geodesies at a specific value of the affine parameter and the 
value of the rate of change of the initial separation at the same value of A, similar to the 
conditions for obtaining particular solutions of second-order, linear, ordinary differential 
equations). 

Now that we have obtained 77 for each case, it is easy to write down the deviation vectors, 
which are given by r? M = rjn^ (where fi = 0,1). Let us now discuss a representative case, say 
E 2 = 4. The range of A (for r > M = 1) is between —00 < A < — 1 and 1 < A < 00. If we 
consider a pair of neighboring geodesies at say A = 1.5 and assume r) is positive, then we are 
forced to choose the solution for rj to be proportional to A 2 . The evolution of rj then suggests 
that as we move further and further to larger positive A, the geodesies spread out (diverge). 
On the other hand, if we assume r) to be negative, then we must choose the \ solution, 
which says that for larger A (which corresponds to larger r as well) , the geodesies converge, 
at least locally. However, as we go to larger A, fj also becomes smaller and approaches zero 
as A — > 00 (r — > 00). Therefore, ultimately the geodesies become parallel to each other. 

It is also useful to note that d 2 r]/d\ 2 , the measure of the gravitational force, is like a 
relative force which acts to change the separation between a pair of neighboring geodesies. 
We can easily substitute the values of A (or A) from the expressions of r(A) in the geodesic 
deviation equations (jUj) , (|2"Bj). and (|28|) to analyze the nature of the force. As an example, 
let us consider the case E 2 = 4. For 77 = A 2 , the relative force is constant throughout and is 
given by: 

Equation (f?TTJ) tells us that the force vanishes as r goes to infinity (recall that the geometry 
is flat in this limit and the geodesies will become parallel to each other), whereas for r — > M, 
it has the value 2/k 3 ^ 2 . A similar analysis can be carried out for the cases with E 2 < 4 and 
E 2 > 4. 

A point worth mentioning here (without proof) is that in the above geometry, the 
geodesies do not seem to converge to a point within a finite value of A in any of the cases 



whereas for 77 = 1/A, we have 
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because of the fact that the Ricci scalar is non- negative (it is zero only when r — > oo). In 
two dimensions, timelike geodesies do not converge to a point (focus) within a finite value 
of the affine parameter unless R < 0.— 

V. CONCLUSIONS 

In summary, the exact solutions of the geodesic deviation equation can provide us with a 
better understanding of the nature of the separation between geodesies. Readers can further 
investigate the deviation vectors in the cases not discussed above in order to improve their 
understanding of geodesic deviation. Trajectories in two dimensions are of the form r(t). 
The nontrivial curvature of the spacetime is responsible for their difference from the usual 
'straight lines' in Minkowski spacetime-the solutions to the geodesic equation obtained in 
the above line element provide examples. If we take pairs of trajectories it is obvious that 
they may converge towards or diverge away from each other-a measure of this effect is 
geodesic deviation. We have illustrated this effect in our example above. 

Our aim has been to provide an example where the geodesic and geodesic deviation 
equations can be exactly solved. It has been shown that for the line element of a two- 
dimensional black hole, exact, reasonably simple solutions do exist. We hope that the 
relevance of the line element, as well as the simplicity of the solutions will attract students 
and teachers to use this example when teaching a first course in general relativity. We 
also mention that, as far as we know, the exact solutions to the geodesic and the geodesic 
deviation equations do not exist in the literature and might be of some interest to researchers 
as well. There are innumerable solutions representing black holes and cosmologies in two- 
dimensional gravity, most of which have appeared largely in the last decade and a half. 
We believe that a carefully chosen section of these line elements can be used as worthwhile 
teaching tools in a first course on general relativity. 
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R. D'Inverno refered to in Ref. |l|. 

The statement made here is based on advanced topics such as the Raychaudhuri equation and 
geodesic focusing theorems (a good, though advanced discussion on these topics is available in 
the text by Wald referred in Ref. 1). In a nutshell, convergence of timelike geodesies (that is, 
focusing) within a finite value of the affine parameter occurs only when the timelike convergence 
condition R^^^ > is obeyed. In two dimensions = ^g^R always and g^v^C = — 1 
(timelike condition). Therefore convergence/focusing occurs within a finite A only if R < 0. For 
some details on this (particularly in two dimensions and in the context of stringy black holes), 
see S. Kar, "Stringy black holes and energy conditions," Phys. Rev. D 55, 4872-4879 (1997). 
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Figure Captions 



Fig. 1: Effective potential (y-axis) as a function of r (x-axis) for three cases E 2 = 1 (Fig. 
la), 4 (Fig. 16), 7 (Fig. lc) respectively; E — 10, k — 1, M — 1 for all three cases. 

Fig. 2: r(A) (y-axis) versus A (x-axis) as in Eqs. 13a (Fig. 2a), 14a (Fig. 26) and 15a 
(Fig. 2c) with A = 3, k = 1, and M = 1. 

Fig. 3: £(A) (y-axis) versus A (x-axis) as in Eqs. 136 (Fig. 3a), 146 (Fig. 36) and 156 
(Fig. 3c) with A = 3, k = 1, and M = 1. 

Fig. 4: r(t) (y-axis) versus t (x-axis) from Eqs. 13a and 136 (Fig. 4a), 14a and 146 (Fig. 
46), 15a and 156 (Fig. 4c) with A = 3, k = 1, and M = 1. 
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Fig. la: E 2 < 4 
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Fig. lb: E 2 = 4 
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Fig.lc: E 2 > 4 

Fig.l: Effective potential (y-axis) as a function of r (x-axis) for three cases E 2 = 1 (Fig. la), 4 
(Fig. 16), 7 (Fig. lc) respectively; Eq = 10, k = 1, M = 1 for all three cases. 



15 



-0.6-0.4-0.2 0.2 0.4 0.6 
lambda 



Fig. 2a: E 2 < 4 




Fig.2: r(A) (y-axis) versus A (x-axis) as in Eqs. 13a (Fig. 2a), 14a (Fig. 26) and 15a (Fig. 2c) 

with A = 3, k = 1, and M = 1. 
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Fig. 3a: E 2 < 4 
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Fig. 36: E 2 = 4 
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Fig. 3c: £ 2 > 4 



Fig.3: t(X) (y-axis) versus A (x-axis) as in Eqs. 136 (Fig. 3a), 146 (Fig. 36) and 156 (Fig 

with A = 3, k = 1, and M = 1. 
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Fig. 4: r{t) (y-axis) versus t (x-axis) from Eqs. 13a and 136 (Fig. 4a), 14a and 146 (Fig. 46), 15a 
and 156 (Fig. 4c) with A = 3, k = 1, and M = 1 in all three cases. 
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